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We present results on the modeling of on- and off-ramps in cellular automata for traffic 
flow, especially the Nagcl-Schreckenberg model. We study two different types of on- 
ramps that cause qualitatively the same effects. In a certain density regime pi ow < 
P < Phigh one observes plateau formation in the fundamental diagram. The plateau 
value depends on the input-rate of cars at the on-ramp. The on-ramp acts as a local 
perturbation that separates the system into two regimes: A regime of free flow and 
another one where only jammed states exist. This phase separation is the reason for the 
plateau formation and implies a behaviour analogous to that of stationary defects. This 
analogy allows to perform very fast simulations of complex traffic networks with a large 
number of on- and off-ramps because one can parametrise on-ramps in an exceedingly 
easy way. 



1. Introduction 



An effective management of transport problems is of great interest in densely popu- 
lated areas. Therefore the performance optimisation of the existing traffic networks 
is an important aspect DLinfra-structural planning §. Recent experimental iHBili] 
and theoretical studies 0'Bu have provided great evidence that the network perfor- 
mance is largely dominated by the capacity of so-called bottlenecks, i.e. parts of 
the roads where the capacity is locally reduced. In realistic traffic systems a large 
variety of possible bottlenecks exists, e.g. road constructions, crossings and lane- 
reductions. If one concentrates on highway networks the reduction of the capacity 
is often due to on- and off-ramps. Therefore the influence of on- and off-ramps has 
been discussed, e.g., on- and off-ramps have been used in studies^ of macroscopic 
models in order to explain the emergence of synchronised traffic 0i3. 
In this work we present a. simu lation study of on^ ari d off-ramps using a discrete 
microscopic traffic model EJJO (for reviews, see lij'Ealia). In order to extract the 
effects of on- and off-ramps we keep the lattice geometry as simple as possible, i.e. 
we study a single-lane system with periodic boundary conditions, with an additional 
on- and off-ramp. Moreover we consider only one type of cars. This choice of the 
system allows for an easy parametrisation and systematic analysis of ramp effects. 
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Nevertheless these simplifications do not reduce the practical relevance of our results 
because the empirical results which have been obtained at multi-lane highways show 



the large influence of the ramps a&B. Moreover the strong coupling of the lanes at 
high densities implies that the generic behaviour is found on all lanes 
The simulation results show that the maximal capacity of the system depends on 
the flow at the ramps as well as on different input strategies. In contrast, different 
output procedures do not affect the results significantly. Starting from this re alis ti c 
scenario we discuss the relation to periodic systems with stationary defects E3'E3, 
i.e. to periodic systems where the mobility of the cars is locally reduced. It turns 
out that the essential features of traffic systems near ramps are already reproduced 
by the system with a stationary defect. Moreover we are able to establish a direct 
relation between on-ramp activity and parametrisation of the defect. Therefore 
it seems to be possible to achieve a considerable reduction of the computational 
complexity of network simulation using this kind of implementation of ramps. 
The influence of on- and off-ramps on the dynamics in microscopic models is not 
yet analysed to our knowledge. The most important point that is investigated here 
is the strategy of the cars to change from the acceleration lane to the driving lane. 
This strategy depends on the individual behaviour of the drivers and might have a 
considerable influence on the traffic dynamics. In this paper we will present a simple 
scenario that is able to quantify the effects of on- and off-ramps qualitatively. It 
will be no restriction to the problem that we will only treat one-lane models. The 
effects caused by the on-ramps in realistic systems can be explained qualitatively 
with one-lane models. Therefore the Nagel-Schreckenberg (NaSch) model serves as 
the basis model for the implementation of on- and off-ramps. Firstly we suggest 
two different types of on-ramps. It is found that it is not necessary to make a 
difference in the modelling of the off-ramps because they play a rather inferior role 
for the traffic dynamics. Afterwards we study the effect of both types on the model 
qualitatively with the help of suitable observables. It will be shown that there exists 
a strict analogy of the effect of on-ramps to that of stationary defects. This analogy 
allows the parametrisation of on-ramps in an exceedingly easy way, so that one can 
perform very fast, but realistic, simulations of complex traffic networks. 
Before we start with our considerations concerning the on- and off-ramps, let us 
remind the reader of the update rules (parallel dynamics) in the NaSch model: 

Rl Acceleration: vi — > mixi{v max , m + 1) 

R2 Braking: Vi — ► min(dj, Vi) 

R3 Randomization: Vi — ► max(0,i>i — I) 

with probability p 

R4 Driving: car i moves Vi cells 



Here Vi and di denotes the velocity of car i and the number of empty cells in front 
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of car i, i.e. the so-called headway, respectively. The maximum velocity and the 
slowdown parameter are denoted as v ma x and p respectively. 

2. Definition of the ramp-types 

In this section we discuss two possible types of implementations of on- and off- 
ramps. The on- and off-ramps are implemented as connected parts of the lattice 
where the vehicles may enter or leave the system. The activity of the ramps is 
characterised by the number of entering (or leaving) cars per unit of time jin (jout)- 
In order to avoid density fluctuations we only added a car to the system if the 
removal of another car at the off-ramp is possible at the same time-step. So, the 
following relation for the rates holds: ji n — jout- Moreover, input and output are 
performed with a constant frequency. Compared to a stochastic in- and output of 
cars this particular choice allows for a better quantitative analysis of the results. We 
have tested that a stochastic in- and output of cars does not lead to a qualitatively 
different behaviour of the system. The chosen length of the on- and off-ramps 
and the distance between them is motivated by the dimensions found on german 
highways. Here we have chosen L ramp = 25 as length of the ramps in units of the 
lattice constant (usually identified with 7.5 meters). The position, i.e. the first cell 
of the on-ramp, is located at x on — 80 and that of the off-ramp at x ff = L — 80 
where L is the system size. Using periodic boundary conditions the distance of 
the on-ramp to the off-ramp is given by d ramp = x on — x Q ff + L. The only and 
essential difference in the implementation of the different types of on-ramps lies in 
the strategy of the cars changing from the acceleration lane to the driving lane. 

Now we discribe two different procedures adding cars to the lattice. 
Type A 

Using this method the lattice will successively be searched in the region of the on- 
ramp (x on to x on + L ramp ) until a vacant cell is found. Then a car will be inserted 
into this cell (see Fig. |l[), even if the cell in front is already occupied by a car. The 
velocity of the car is set to v max . Dependent on the global density p this can lead 
to a strong perturbation of the system (see Section ). 

This extremely restrictive method described above implies a very inconsiderate be- 
haviour of the drivers because in the process of changing from the acceleration lane 
to the driving lane no safety margin will be held. However, due to the unbounded 
deceleration capabilities of the cars in the NaSch model this will not lead to acci- 
dents. 
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Type B 

In this method the cell that will be occupied by a car is selected in a stochastic way. 
A vacant cell will be randomly chosen in the region of the on-ramp. Afterwards a 
car will be inserted into this cell like in type A. The measurements in the simulations 
will show that there is no qualitative difference of this type to type A. Type B will 
also lead to a strong perturbation of the system. 

By the method described above one gets a more realistic description of traffic flow 
than in type A. Here it will be used to simulate both inconsiderate and cautious 
behaviour of the drivers. 



On-ramp 



Inserted car 



Figure 1: Schematical representation of the on-ramp (Type A). 

The off-ramps work for the corresponding types of on-ramps in the same way: one 
goes successively through the cells of the lattice in the region of the off-ramps until 
an occupied cell is found. Then the car will be removed from this cell. A schematic 
diagram of this procedure is shown in Fig.^. 

Off-ramp 



Removed car 

Figure 2: Schematical representation of the off-ramp. 



3. Effects of the ramps 

In this section we discuss the effects of the ramps. Using computer simulations 
fundamental diagrams dependent on the in/output-rates and density profiles de- 
pendent on the global density of the system are generated. Then, by analysing 
these results, we investigate the relation between on-ramps and stationary defects. 
This will be done in the last subsection. 
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3.1. Fundamental diagrams and density profiles 

In Fig. ^ the fundamental diagrams of the NaSch model for both types of on-ramps 
are shown. The input-rate is chosen to be j in — j out — i. This choice guarantees 
the conservation of the global density p. For comparison the fundamental diagram 
without ramps is also shown as solid line. 
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Figure 3: Fundamental diagrams of both ramp types at the same input-rate ji n = g. 
The model parameters are given by L = 3000, v max = 5,p = 0. 

It is clearly seen that a density regime pi ow < p < Phigh exists where the flow J(p) is 
independent of the density. This so-called plateau value of the flow is lower than the 
corresponding flow of the model without ramps. An increase of the input-rate jin 
leads to a decrease of the plateau value If one compares the different types of input 
strategies it is evident that the plateau value of type B is lower than that of type 
A. This can be explained by the probability that the cell in front of the inserted 
car is already occupied. For type B this probability is smaller than for type A. 
Nevertheless there is no difference between the two types in a qualitative sense. 
That is the reason why we restrict ourselves to type A in the further treatment. 
With this type the effect can be seen in the simplest way and so it is easier to 
describe. 

Obviously we can distinguish three different phases depending on the global density. 
In the high and low density phases the average flow J(p) of the perturbated system 
takes the same value as in the system without ramps. For intermediate densities 
Plow < P < phigh the flow is constant and limited by the capacity of the ramps. 
For a better understanding of the behaviour of the average flow it will be helpful to 
look at the density profile (see Fig. ||). In the high and low density phase one only 
observes local deviations from a constant density profile. For intermediate densities 
the system is separated into macroscopic high {phigh) and low (pi ow ) density regions. 
In the region of the ramps the density p ram p is additionally higher than in the high 



NaSch 
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density region. Thus the ramps act like a blockage in the system that decreases the 
flow locally. Varying the global density within the phase-separated state the bulk 
densities in the high and low density region remain constant, only the length of the 
regions changes (see Fig. ||). Also the local density at the ramps (p r amp) remains 
constant. Finally this leads to a constant flow in the segregated phase. 





p=0.1 - 




p 0.3 




p=0.6 


















L 


_J 





500 1000 1500 2000 2500 3000 



Figure 4: Density profiles in the three phases. In the high (p = 0.60) and low 
(p = 0.10) density phases only local inhomogeneities occur near the ramps (x on = 
80, %off = L — 80). For intermediate densities (p — 0.30) one observes phase 
separation. The model parameters are the same we used in Fig. |J 
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Figure 5: Density profiles for two different densities within the phase-separated 
state. Pramp, Plow, Phigh remain constant, only the length of the regions changes. 

If one assumes that the cars can move freely in the low density regime and that 
the flow in the high density region depends linearly on the density (which holds for 
small braking parameters p) we get the following relation: 

Plow = (1 - Phigh) • (3-1) 

Vmax P 

Unfortunately no expressions for pi ow and Phigh depending on the model parameters 
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only can be given. The reason is that no exact analytical description for the unpcr- 
turbated system with v max > 1 has been found so far. Furthermore an analytical 
treatment of the ramps is very complicated because of the strong interactions within 
these regions. 



3.2. Analogy to stationary defects 

In the previous subsection we have shown that the on-ramps in the system act like 
a blockage. This blockage leads to a decrease of the local flow at the ramp. To 
clarify the analogy of on-ramps to stationary defects we give a short summary of 
the most important results on stationary defects like construction sites EUll. 




Figure 6: Schematical representation of a single defect in the NaSch model with 
Vmax = 1- In this special case the average velocity of the cars is approximately 
v = 1 — p for p <C i. At the defect site the average velocity is v = 1 — pd- 

Stationary defects can be implemented in two different ways. First, one can define a 
certain range on the lattice where th e sl owdown parameter p is increased compared 
to that of the residual lattice sites_p. Second, one can decrease the maximum 
velocity v max in this certain range r°H 21 l. In our investigations we restrict ourselves 
to the first method. 

Fig. H shows a schematical representation of a single defect in the NaSch model 
with v m ax — 1 and slowdown parameter p. The slowdown parameter at the defect 
is denoted as Pd- 

In Fig. ^ the fundamental diagram of the NaSch model with stationary defects is 
represented. Again we can distinguish three different phases. In the high and low 
density regime the flow J(p) is identical to that of the model without defects. For 
intermediate densities pi ow < p < Phigh the flow is constant in analogy to the model 
with ramps. In this regime J(p) is limited by the capacity of the defect sites. 
Once again the phenomenon of the plateau formation can be understood by the 
analysis of the density profiles (see Fig. ||) . The argument is the same as in subsec- 
tion 3.1. 

For the case v max = 1 of the NaSch model with a single defect site good approxi- 
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Figure 7: Fundamental diagram of the NaSch model with stationary defects. The 
model parameters are given by L = 3000, p = 0,v max = 5. For comparison the 
diagram of the model without defects is shown as solid line. 



mations for the plateau value of the flow as well as for the densities pi ow and phigh 
can be made cB. An essential reason for this is the particle-hole symmetry in this 
special case. Because of this symmetry following relation holds: 

Plow Phigh 1. (3.2) 

With the assumption El that the system is separated into two regions of constant 
density and that the flow Jd at the defect is equal to the flow Jbuik in the bulk one 
can establish following equation in the mean field theory: 

Jd = QdPhigh (1 — Plow) = Jbulk = QPhigh (1 — Phigh) (3.3) 



where q = 1 — p and qa = 1 — Pd- With equation fl3.2| ) we finally get the expressions 
for the densities: 

_ q Id /„ .s 

Phigh — — ; > Plow — — ; • 

q + qd q + qd 



From this and the exact current (see e.g. tBEn) follows the plateau value of the flow: 



Jp 



Ul--^^/(q + qdf-M 2 qd). (3.5) 
2 \ q + qd * ) 



In principle one can treat the defects for v max > 1 in the same way. When we make 
the assumption that for the two branches in the fundamental diagram the following 
approximations hold: 



J free = P (Vmax ~ P) , Jjam = (1 — p) (1 — /)) , 



(3.6) 
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Figure 8: Density profiles of the NaSch model with stationary defects for densities 
in the three phases. In the high (p = 0.60) and low (p = 0.10) density phases only 
local inhomogeneities occur near the defect sites (L — v max < Xd < L). For densities 
within the plateau regime one observes phase separation. 



one gets the following expressions for the densities in the phase-separated state: 

Qd 



Pi 



OW 



Vraax ~ P + 9d ' 



(v m ax - p)(q - qd) + qqd ,„ 

Phigh = 7 r ; • •>■' 

{Vmax - P)q + qqd 
The plateau value of the flow results from this as: 

{ v max P)qd 

Vmax - P + qd 

In spite of the crude estimate for Jj am this approximation is in good agreement 
with the simulation data for small p (see Fig. ^) . It has not been checked so far for 



which range of values for pd equation (3.8) can reproduce the simulation results 



We have shown that there is no qualitative difference between the effect of on- 
and off-ramps and that of stationary defects. In both cases one observes plateau 
formation in the fundamental diagram as well as phase separation in the system. 
The only difference lies in the nature of the blockage dividing the system into two 
macroscopic regions. In the case of the ramps it is the local increase of the density 
that decreases the flow locally. In the model with defects the increased slowdown 
parameter leads to a local decrease of the flow. 

In order to use the equivalence of ramps and defects in realistic simulations of traffic 
flow one has to determine the relation between the input-rates jin and the slowdown 
parameter pd at the defect sites. Fig. || shows the values of pd and jin leading to 
the same fundamental diagrams. It can be seen that a non-trivial relation between 
the two parameters exists. We regard this Fig. || as our main result for future 
applications. 
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Figure 9: Relation between the input-rates ji n and the slowdown parameter pd at 
the defect sites. 



4. Summary and Conclusions 

In this paper we have presented a simple scenario that can quantify the effects of 
ramps on the traffic dynamics in microscopic models. The model used here is the 
well known Nagel-Schreckenberg model, a cellular automaton for traffic flow. We 
have introduced two different types of on-ramps. One type that implies a very 
inconsiderate behaviour of the drivers and another type that simulates both incon- 
siderate and cautious behaviour. We have shown that there is no difference between 
these two types in a qualitative sense. With both types one observes effects very 
similar to those of stationary effects: in both cases one observes plateau formation 
in the fundamental diagram as well as phase separation in the system. The only 
difference lies in the mechanism dividing the system into two macroscopic regions. 
This analogy is is very important for efficient modelling of complex networks El with 
a multitude of ramps. 

Usually in simulations of large networks one uses a more realistic approach. In El 
ramps are fed by an external road. Cars try to change to the highway according to 
the usual lane changing rules. If they could not change the cars wait at the end of 
the ramp. 

In our considerations we have also investigated in the influence pf ramps on different 
mutations of the NaSch model like the Slow To Start model oEj and a model with 
an anticipation rule. The effects caused by the ramps are qualitatively the same 
as in the NaSch model, but in the model with anticipation one has to use much 
higher input-rates j m to observe phase separation and plateau formation. This 
model reacts more robust on the perturbations caused by the cars changing from 
the acceleration lane to the driving lane cl. 

The effects of ramps have been studied in connection with synchronised traffic 
recently Bi. For the hydrodynamic models it has been found that the localised 
perturbation through the ramp flow can induce transitions between different traffic 
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states. 

In □ it has been shown using computer simulations and data from measurements on 
real traffic that on-ramps can induce a first-order non-equilibrium phase transition 
between free flow and a congested phase. The transition is induced by the interplay 
of density waves caused by the on-ramps and shock waves moving on the highway. 
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